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Abstract 

We present lightcurve observations and multiband photometry for 
107P/ Wilson-Harrington using five small- and medium-sized telescopes. The 
lightcurve has shown a periodicity of 0.2979 day (7.15 hour) and 0.0993 day 
(2.38 hour), which has a commensurability of 3:1. The physical properties 
of the lightcurve indicate two models: (1) 107P/Wilson-Harrington is a 
tumbling object with a sidereal rotation period of 0.2979 day and a precession 
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period of 0.0993 day. The shape has a long axis mode (LAM) of Li:L2:L^ 
= 1.0:1.0:1.6. The direction of the total rotational angular momentum is 
around A = 310°, /3 = -10°, or A = 132°, /S = -17°. The nutation angle is 
approximately constant at 65°. (2) 107P/Wilson-Harrington is not a tumbler. 
The sidereal rotation period is 0.2979 day. The shape is nearly spherical but 
slightly hexagonal with a short axis mode (SAM) of Li:L2:L^ = 1.5:1.5:1.0. 
The pole orientation is around A = 330°, /3 = —27°. In addition, the model 
includes the possibility of binary hosting. For both models, the sense of 
rotation is retrograde. Furthermore, multiband photometry indicates that 
the taxonomy class of 107P/Wilson-Harrington is C-type. No clear rotational 
color variations are confirmed on the surface. 
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1. Introduction 



Asteroids and comets are primordial bodies that formed in the earliest 
stage of the solar system. Their rotational states, shapes, and material 
reflect the collisions, disruptions, and chemical processes since then to 
the present. Some small solar system bodies exhibit behavior such 
as that shown by both comets and asteroids (so-called, comet-asteroid 
transition objects). As an example, near-earth object (NEO) (3200) 
Phaethon shows signs of past cometar y activity b ecaus e it is thought 
to be associated with the Geminid fl Gust af son! . Il989[ ). Dynamical 
numerical simulations and spectral observations for (3200) Phaethon 
support (2) Pallas, which is outer m ain belt as t eroid s, is the most 



likely parent body of (3200) Phaethon ( IClark et al. 



20101 



de Leon et al. 
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20101). Meanwhile, objects that display cometary activities in the 



main-belt asteroid (MBA) region have recently been discovered . They are 



classified as main-belt c omets (MBCs) f|Hsieh and Jewit 



are 133P/Elst-P izzaro flElst et al.l. Il996h . P/2005 Ul fiRead et al.l. l2005h 



176P/LINEAR flHsieh et al 



2006): the MBCs 



201 If). P/2008 Rl (Garrad) flJewitt et al. 



2009|), P/2010 A2 fL 



Sagral (jMarsden et al. 



Jewitt et al. 



N 



EAR ) flBirtwhistle et al.l . hom. P/2010 R2 fLa 



20101), and (596) Scheila (IBodewits et al. 



2011 



20111 ). One possible activat ion mechanis r n for MBCs is impacts 



with s mall (e.g., meter-si z ed) objects (iTothl . 



2008 




Jewitt et al. 


Jewitt et al.l. 


2011) 



2010 



Snodgrass et al.l . 



2000 



20101: 



Diaz and Gil-Hutton . 



Bodewits et al. 



2011 



20111 ). The other activation mechanisms are rotational- fissions 



due to the spin-up by Y arkovsky-Q 'Keefe-Radzievsk ii-Paddack 



effects ( Jewitt et al. 



(YORP) 



20101 ). and thermal infiuences (I Jewitt et al. 



20091). 



Interesting properties of MBCs are their dynamical origin and possible 
fun ction as reserv oi rs for water-ice and organics. A numerical integration 



by 



Haghighipourl f|2009[ ) states that the origin of 133P/Elst-Pizarro, 



176P/LINEAR, and P/2005 Ul (Read) is concordant with the Themis family 
of asteroids. Compared with all asteroids, the Themis family of asteroids 
includes B-type asteroids at a relatively high population rate. Some B-type 



asteroids in the T h emis f a mily seem to 



(lYang and Jewitt 



2010 



Clark et al. 



lave experienced aqueous alterations 



20101 ). (3200) Phaethon is also a 



■ 3-type asteroid and shows the existence of aqueous alteration materials 



Licandro et al 



20071 ). In addition, water-ice and o rganics are detected on 



the surface of (24) Themis (IRivkin and Emeryl . |2010[) and (65) Cybele , which 
orbits along the outer edge of the main belt (ILicandro et al.l . l201ll ). The 
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study of comet-asteroid transition objects provides keys to the dynamical 
origin and evolution of NEOs, the mutual collisions of small solar system 
bodies, the material differences between asteroids and comets, and the origin 
of Earth's water. 

This study's purpose is to obtain the rotational states, shape model, 
and rotational color variations for 107P/Wilson-Harrington (also know 
as (4015) Wilson-Harrington; hereafter 107P), which is a representative 
comet-asteroid transition object. 107P was discovered accompanied by 



a fai nt cometary tail at Palomar Observatory in 1949 (IFernandez et al 



19971 ). The object, however, could not be tracked because of insufficient 



observations to determine an accurate orbit. Later, a near-earth asteroid 
1979VA (= 4015) was discovered. Subsequent observations identified asteroid 
(4015) 1979VA and 107P as the same object. Despite a devoted search, no 



cometary activity has been detected since the initia . 



( IChamberlin et al. 



1996 



Lowry and Weissmanl . 



2003 



observation o: 



Ishiguro et al. 



107P 



201l|). 



107P is an Apollo asteroid whose orbital parameters are a = 2.639 AU, e = 
0.624, i = 2.7 85°, and the Tisser and parameters (Tj) = 3.08. A numerical 



simulation by 



Bottke et al 



( 120021 ) mentions that there is a 4 % chance that 



107P has a JFC origin and a 65% chance it has an origin in t he outer 



main-belt region. Taxonomically, it is categorized as a CF-type (ITholenl . 



19891 ). The refie ctance spectrum in the region 3800-6200 A is similar to 



(3200) Phaethon (IChamberlin et al.l . 119961 ). The thermal properties of 107P 
have been investigated by mid -infrared photometry with NASA's Spitzer 



Telescope (ILicandro et al 



20091 ). These observations show that the beaming 



parameter, the diameter, and the albedo are rj = 1.39 ± 0.26, D = 3.46 ± 0.32 



4 



km, and p„ = 0.059 ± 0.011, respectively. The rotational period of 107P has 



been r eport ed to be 3. 5 56 ho ur and 6.10 ± 0.05 hour b y 



f ll983h and 



Osip et al. 



Osip et al. 



Harris and Young 



fll995h ascri 



Des the 



Harris and Young 



fll983[ ) 



(Il995[ ). respectively, 
difference of the two reports to the noisy data of 
because of the weather conditions. The few days' observation in both reports, 
however, is not enough to determine the correct rotational period. Longer 
observations are required to derive the correct rotational period and other 
physical properties. 

We hypothesize that 107P migrates to the NEO region from the outer 
main-belt region inhabited by six of seven known MBCs, and impacts with 
small objects could eject dust and/or expose sub-surface ice that then trigger 
107P's cometary activity. Post-MBC, 107P is capable of becoming host to 
water-ice, organics, and aqueous alteration materials. In this hypothesis, 
the impacts' influence would be apparent in the rotational states and/or the 
surface color variations. 

We had an opportunity to observe 107P from August 2009 to March 2010. 
Our long observation campaigns enable us to derive the rotational states, 
shape model, and rotational color variations. Furthermore, the orbit of 107P 
makes it accessible by spacecraft. A more advanced sample return mission 
from a D-type asteroid or an ast eroid-comet transition o biect is envisioned in 
Japan. One candidate is 107P (lYoshikawa et al.l . |2008| ). Clarification of the 
physical properties of 107P is important to the design of the future mission. 



If we are able to obtain 107 
revise the physical model of 



s physical propertie s, the data will be useful to 



Licandro et al 



(|2009[ ) , similar to the Hayabusa-2 



target 162173 (1999JU3) whose physical model was reconstructed by both 
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thermal observations and the hghtcurve (iMiiller et al.l . 1201 If ). This paper is 
organized as follows. In Section 2, we describe the observations made and 
the data reduction. In Section 3, we mention the rotational states and shape 
model of 107P. In Section 4, we focus on the possibility of tumbling motion 
and the existence of a binary. Finally, we summarize the physical model of 
107P and discuss the feasibility of a sample return mission. 

2. Observations and Data Reduction 
2.1. Observations 

We conducted the photometric observation campaigns of 107P with 
five small- and medium-sized telescopes. The observational circumstances 
and the states of 107P are listed in Table 1 and Table 2, respectively. 
All telescopes were operated with the non-sidereal tracking mode. The 
longest-term observation of this campaign was carried out using the 1.0 m f/3 
telescope at the Bisei Spaceguard Center (BSGQH) from September 6, 2009 
to March 11, 2010. The detector consisted of four CCD chips with 4096 x 
2048 pixels. We used one CCD chip to obtain as many images as possible by 
shortening the processing time. The field of view (FOV) for one CCD chip 
is 1.14° X 0.57° with a pixel resolution of 1.0". The exposure time varied 
from 30 s to 600 s according to the observational situations. Individual 
images were taken with a commercially available short-pass (long-wavecut) 
filter, the effective wavelength of which ranged from 490 nm to 910 nm. 
We denote the filter as W in Table 1. In order to investigate rotational 



^BSGC is administrated by the Japan Space Forum. 
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color variations, multiband photometry was conducted using a Sloan Digital 
Sky Survey (SDSS) g' , r', i' , z' filter on December 17, 2009. One set of 
observations was made using three consecutive images for each filter. The 
filters were changed in the following sequence: three g' images — >■ three r' 
images — )■ three i' images three z' images. We repeated this sequence four 
times. 

The second-longest-term observation used the 0.5 m f/6.5 Multicolor 
maging Tele s cope for Survey and Monstrous Explosions (MITSuME) 



Kotani et al. 



20051 ) at Okayama Astrophysical Observatory (OAO) from 



November 7, 2009 to December 21, 2009. The telescope is capable 
of obtaining a three-color (SDSS g' , Johnson-Cousins Rc and Jc) image 
simultaneously. The detector is 1024 x 1024 pixels CCD with FOV of 26' 
X 26' (1.52"/pixel). The images were taken with an exposure time of 120 s. 
To search for the rotational color variation, we used the data of December 
17 because they could be compared with the observations of the BSGC and 
the photometric precision of the other day's data was not sufficient to detect 
the color variation. 

The third observation was carried out using a 1.05 m f/3.1 Schmidt 
telescope with 2048 x 2048 pixels CCD at Kiso Observatory on August 17, 
19, and 20 and December 12, 2009. This instrument provides a FOV of 50' x 
50' (1.46"/pixel). The images were obtained using a Kron-Cousins Rc filter 
with an exposure time of 120-300 s. 

The fourth observation was made using the Lulin One-meter Telescope 
(LOT) Jfluang et al.l . 120051 ) in Taiwan on December 7-10, 2009. The CCD 
consists of a 1340 x 1300 array, and the FOV covers the area of 11.5' x 11.2'. 
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The pixel resolution and f-number are 0.51"/pixel and 8, respectively. The 
images were obtained using a Johnson-Cousins Rc filter with an exposure 
time of 90 s. 

The last, observation was made using the University of Hawaii 2.24 m 
f/10 telescope (UH88) on December 19, 2009, with a 2048 x 2048 pixels 
CCD. The FOV of the instrument is 7.5' x 7.5' with a pixel resolution of 
0.44". Almost all images were obtained using a Kron-Cousins Rc filter with 
an exposure time of 60 s. 



Table 1: Observation states. 



Observatory 


Year/Mon/Day 


Exp time(s) 


Filter 


BSGC 


2009/12/17 


300 


g', r', i', z' 


BSGCt 


2009/09/6,7,9,10,15,16,19, 2009/10/8,10,28 


30-180 


W 


BSGC 


2009/11/3,5-7,11,14, 2009/12/5,7-9,19,22 


60-300 


W 


BSGC 


2010/01/3,6-8,14-18,22,23 


180-300 


W 


BSGCt 


2010/02/3-5,7,9,16,18,19, 2010/03/11 


240-600 


W 


OAO 


2009/11/7,14,15,18-21,23, 


120 


g', Rc, Ic 


OAO 


2009/12/1,2,6,7,14,16-21 


120 


g', Rc, Ic 


KISOt 


2009/08/17,19,20, 2009/12/12 


120-300 


Rc 


LOT 


2009/12/7-10 


90 


Rc 


UH88 


2009/12/19 


60 





t A sufficient number of data was not obtained from August to October because the altitude of 107P fell 
below 25° about 30 minutes from the observation start and 107P overlapped stars of the galactic plane. 
The photometric precision was insufficient after February 7, 2010. We did not use these data for the 
estimation of rotational periods and shape models. These data were utilized for the trend confirmation of 
lightcurves that were obtained from the other day's data and the monitoring of cometary activity. 
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Table 2: States of 107P in each month. 



Ycar/Mon/Day 


A [AU]'' 


R [AU]'' 


a 


o jc 


Sky motion ["/min] 


m 


d 


2009/08/17-20 


0.687-0.684 


1.309-1.286 


49.8- 


-51.4 


0.45-0.55 


17.7- 


17.7 


2009/09/6-19 


0.653-0.612 


1.162-1.084 


59.9- 


-66.0 


1.09-1.49 


17.9- 


17.9 


2009/10/8-28 


0.529-0.434 


1.008-0.995 


73.8- 


-77.2 


2.18-3.26 


17.3- 


17.5 


2009/11/3-23 


0.410-0.382 


1.006-1.083 


76.3- 


-65.4 


3.65-4.52 


17.1- 


16.7 


2009/12/1-22 


0.401-0.543 


1.130-1.277 


59.2- 


-46.3 


4.31-2.96 


16.6- 


17.0 


2010/01/3-23 


0.670-0.936 


1.373-1.539 


42.0- 


-37.8 


2.32-1.70 


17.4- 


-18.4 


2010/02/3-19 


1.105-1.372 


1.632-1.767 


36.1- 


-33.8 


1.51-1.34 


18.8- 


19.7 


2010/03/11 


1.732 


1.932 


30.8 


1.22 


20.0 



^ Object to observer distance. 
^ Heliocentric distance. 
Phase angle (Sun-107P-obscrver). 

Apparent magnitude. This value is estimated using UCAC 2 catalog stars that arc taken in the same 
field with 107P. 

2.2. Data reduction 

All images were bias and flat-field corrected. When using the data of 
OAO for the derivation of lightcurve, we stacked four images to compensate 
for the poor flux. All observation times were corrected using the light travel 
time from 107P to the Earth. By using the IRAF/APPHOTu package, we 
measured the raw magnitude of 107P and from three to seven reference stars 
that were bright enough compared with 107P. We set aperture radius to 1.5 x 
FWHM for 107P and reference stars images, respectively. Since the reference 
star images are slightly elongated by the non-sidereal tracking, the aperture 



^IRAF is distributed by the National Optical Astronomy Observatory, which is 
operated by the Association of Universities for Research in Astronomy (AURA) under 
cooperative agreement with the National Science Foundation. 
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radius is larger than that of 107P. We cahbrated the magnitude fluctuations 
due to the change of atmospheric conditions as follows, 



F^{t) = F^{t)-F^{t). 



(1) 



Here, F^{t) is the lightcurve by rotation of 107P in i-th observation day, 
F^{t) is the raw magnitude of 107P, F^{t) is the averaged raw magnitude of 
reference stars and represents the change of atmospheric conditions, and t is 
the observational time. Next, we define the averaged magnitude of F^{t) in 
each night as the normalized (zero) magnitude. The lightcurve by rotation 
of 107P can be rewritten as 



where F^ is the averaged magnitude of Flit). Since the averaged magnitude 
is normalized to zero magnitude for all nights, we can connect the different 
night's lightcurve with little regard for the difference of absolute magnitude. 
In addition, the difference of reference stars each night does not affect the 
periodicity of lightcurve. The problem of this procedure could include the 
offset between different nights when the short observation time per day poses 
the detection of a specific peak (bottom) in the lightcurve. However, the 
peaks and bottoms in the lightcurve have been detected evenly (See Fig. 3) 
because the observation time per day is long enough from November 2009 to 
February 2010 when the data are utilized for the analysis. Thus, the offset is 
negligible. Furthermore, the apparent magnitude change of 107P is gradual 
up to 1.0 magnitude per month (Table 2). The change does not act on the 
derivation of rotational period that is expected to be from three to seven 



Fl^^it) = Flit) - Fl, 



(2) 
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hours according to past reports. 

In contrast to the relative photometry of the hghtcurve, more photometric 
precision is required to detect the rotational color variation by multiband 
photometry. In order to improve the photometric precision, we averaged 
three consecutive images of 107P for the BSGC's data and 14-16 consecutive 
images of 107P for the OAO's data. We also measured the flux of ten 
standard stars from SDSS data Release 7 ( Abazajian et al.l . 2009 ). whose 
stars were imaged simultaneously in the same frame as 107P (Table 3). These 
objects have magnitudes of about 14 mag to 16 mag in the r'-band and 
classification code 1 (= primary), quality flag 3 (= good), and object class 
6 (= star). We evaluated atmospheric extinction coefficients and conversion 
factors to standardize the SDSS system for each filter. The atmospheric 
extinction coefficient was calculated by the magnitude variations of the 
standard stars for the change in airmass. Extra-atmospheric instrumental 
magnitudes of both 107P and the standard stars were derived using the 
obtained atmospheric extinction coefficient. The conversion factors were 
estimated by comparing the extra-atmospheric instrumental magnitudes with 
the magnitude of standard stars. In BSGC's observation, the multi-color 
images were not obtained simultaneously. The brightness of 107P by the 
rotation changes inevitably during the filter switch. We defined the time of 
the third r' images in each sequence as a standard time, and then calibrated 
the amount of brightness change for the standard time. The amount of 
brightness change was estimated by the fitting curve of the lightcurve. In 
order to compare the OAO's data with BSGC's, the Rc and Ic magnitudes 
obtained at OAO were converted to r' and i' magnitudes using the conversion 
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equations proposed by iJordi et al.l (120061 ) . 

Table 3: Standard stars in SDSS-7. 



Ra[° 



Dcc[°l 



12.051081 


8.615352 


15.891 ± 0.003 


15.279 ± 0.003 


15.057 ± 0.003 


14.913 ± 0.005 


12.246608 


8.604861 


15.209 ± 0.003 


14.756 ± 0.003 


14.598 ± 0.003 


14.525 ± 0.004 


12.073952 


8.515561 


14.536 ± 0.003 


14.215 ± 0.003 


14.115 ± 0.003 


14.074 ± 0.004 


12.011267 


8.489718 


15.766 ± 0.004 


14.891 ± 0.004 


14.587 ± 0.003 


14.444 ± 0.004 


12.275938 


8.545386 


16.158 ± 0.003 


15.566 ± 0.003 


15.349 ± 0.004 


15.230 ± 0.005 


11.898784 


8.539119 


14.947 ± 0.003 


14.554 ± 0.003 


14.418 ± 0.003 


14.366 ± 0.004 


12.365509 


8.569521 


14.587 ± 0.003 


14.037 ± 0.003 


13.839 ± 0.003 


13.713 ± 0.003 


12.374946 


8.538940 


16.005 ± 0.003 


15.555 ± 0.003 


15.396 ± 0.004 


15.307 ± 0.005 


11.796300 


8.569376 


15.308 ± 0.003 


14.641 ± 0.003 


14.413 ± 0.003 


14.287 ± 0.004 


12.371825 


8.478789 


15.763 ± 0.004 


15.298 ± 0.004 


15.156 ± 0.004 


15.099 ± 0.005 



3. Results 

3.1. Rotational states 

Since the "Standard Feature (SF)" that is, the flux peaks and/or bottoms 
in hghtcurves, shifts along the phase of lightcurves due to changes in 
the geometric relationship between the Earth, 107P, and the Sun during 
the long-term observations, the estimation of the sidereal rotation period 
for 107P requires a short-term observation within a few weeks. We use 
the data from December 7 to 22, when the phase-shift is small. The 
photometric precision of 107P and the observational implementation time 
per day are enough to make clear the sidereal rotational period during the 
term. Assuming double-peak ligh tcurves , a period ana l ysis i s carried out 



1976 



Scargld . Il982l ). The power 



with a Lomb-Scargle periodgram (iLombl . 
spectrum from the period analysis shows four period candidates of 0.0993 
day, 0.2294 day, 0.2591 day and 0.2979 day (Fig. 1). A typical error of 0.0002 
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.A/\a/v.Aa/^ 



0. 



0.: 



0.4 0.6 
Period[day] 



0.8 



Fig. 1: Power spectrum for the sidereal rotation period of 107P by assuming the 
double-peak lightcurve. The calculation is carried out on data obtained from December 7 



day corresponds to ± 0.005 hour. Though the most significant candidate is 
0.2591 day, we conclude that 0.2979 day 7.15 hour) is the sidereal rotation 
period of 107P for the following reasons. First, the amount of the amplitude 
in the folded lightcurve with 0.2591 day is not stable in the same phase. That 
is to say, the different amplitudes overlap on a specific phase. For example, 
the lightcurve peaks and bottoms overlap around the phase of 0.2-0.4 in the 
folded lightcurve with 0.2591 day (Fig. 2: Top). In the case of the folded 
lightcurve with 0.2979 day, the same amplitudes appear periodically (Fig. 
2: Bottom. A few lightcurves each night are also shown in Fig. 3). The 
periods of 0.2591 day and 0.2979 day correspond approximately to 3.86 and 
3.36 cycles per day, respectively. The difference is just 0.5 cycles per day. 
Lightcurves mainly represent the light scattering cross section of objects. 
When we assume that an object is a symmetric ellipsoidal body, almost the 
same cross section appears in every half rotation. Therefore, it is difficult 



to 22. 
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Fig. 2: Lightcurve of 107P. (Top) The liglitcurve is folded with 0.2591 day. The peak and 
bottom of the lightcuve overlap around the phase between 0.2 and 0.4. (Bottom) The 
lightcurve is folded with 0.2979 day. The same amplitudes appear periodically. 
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Phase 

0.4 0.6 0.8 1 1.2 



-0.3 - 




0.3 - 

172.85 172.9 172.95 173 173.05 173.1 173.15 
JD-2455000[day] 
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-0.3 - 
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0.2 - 
0.3 - 

173.85 173.9 173.95 174 174.05 174.1 174.15 
JD-2455000[day] 



Fig. 3: (Top) Lightcurve in December 7, 2009. (Bottom) Lightcurve in December 8, 2009. 
The phase corresponding to Fig. 2 is added to the top scale of the figures. 
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to distinguish the difference of the half rotation using the short observation 
time, which is comparable with the sidereal rotation period. We call the 
indistinctive period a pseudo-period. The periods of 0.2591 day and 0.2294 
day (= 4.36 cycles per day) are the pseudo-period of 0.2979 day. Second, 
the period of 0.2979 day is able to explain the previous reports about the 
sidereal rotation per iod of 107P. Since the p eriod of 0.2979 day is around 
twice the period of [Harris and Young] (119831 ) (0.1482 day ~ 3.556 hour), 
their data show enough periodicity in 0.2979 day. Needless to say, assuming 
the lightcurve of 107P has a triple-peak, the period of 0.1490 day (~ 3.58 
hour) is also a candidate for the sidereal rotation period. However, the 
possibility of 0.14 9 day is eliminated by th e inconsistency with the data of 



Osip et al 



(119951 ). iHarris and Young! (119831 ) would recognize their lightcurve 
as the typical double-peak with the period of 3.556 hour., because the third 
a mplitude of fl ux in their lightcurve was not detected. Moreover, the period 



of lOsip et al.l (jl995l ) (0.2542 day = 6.1 ± 0.05 hour) is approximately the 
same as 0.2591 day (~ 6.22 hour). Our data set also has a sufficiently high 
significance level around the period of 0.2542 day. As we mentioned above, 
however, the perio d of around 0.25 91 day is a pseudo-period. Since the 

(119951 ) was only two days, the demarcation of 



observation term of 



Osip et al 



a pseudo-period would have been difficult. Third, we focus on the lightcurve 
of 0.2979 day as having an unusual six peaks. The period of 0.0993 day is just 
one third that of 0.2979 day. If 107P has a typical double-peak lightcurve, the 
period of 0.0993 day is the sidereal rotation period. However, the amplitudes 
overlap the different three peaks and bottoms in the folded lightcurve with 
a period of 0.0993 day. Thus, we exclude the period of 0.0993 day as the 
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sidereal rotation period. 

On the other hand, the period of 0.0993 day may be the precession 
period. If an object has tumbhng motions, the hghtcurve is dominated by 
two periods: one, P^, for the rotation about the extremal axis of the object 
as an inertia ellipsoid, and the other, P^, for the precession about the total 
rotational angular momentum vector. When frequencies are defined as 2/^ 
= Pr^, the lightcurve periodicity of tumbling ob jects appears 



P^' and 



at fr equencies that are a linear combination of and (iKaasalainenl . 
2OOII ). When we assume that P^ is 0.2979 day and is 0.0993 day, the 
frequency 4/^ = % fd> — fih) = 6.713 day~^ ap proximately corresponds to the 



inverse of period of 



Harris and Youngl (jl983[ ) and one half of our rotational 



period of 0.2979 day. The existence of periodicity of the linear combination 
of two periods shows circumstantial evidence for tumbling. We make the 
shape model of 107P in the following subsection and discuss the feasibility 
of tumbling motion in Section 4.1. 

3.2. Direction of total rotational angular momentum and shape model 

Above, we suggested the possibility of tumbling motion. If 107P is a 
tumbler, the pole orientation does not accord with the direction of total 
rotational angular momentum and is not stable. What we can obtain 
is not the pole orientation but the direction of total rotational angular 



the "epoch method" ( 


Maenusson 


1986 


) or the "lightcurve 


(Kaasalainen and Torppa. 


2001 




Kaasalainen et al.. 


2001. 



20021). The "amp litude method" is also proposed as an alternative method 
( lMagnussonlll986[ ) . However, we cannot adopt the amplitude method because 
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of the small amplitude change during the observational term. The epoch 
method determines the direction of total rotational angular momentum by 
minimizing the phase-shift of the SF in lightcurves. We select a lightcurve 
peak around the phase of 0.01 at the bottom of Fig. 2 as the SF because the 
peak is better observed than any other feature. Moreover, we use the data 
obtained from November 5, 2009 to February 5, 2010. The identification of 
the lightcurve peak is difficult from the data of another term because of the 
photometric error. The phase-shift can be written as 



where Tq is the time at the first SF, Ti is the time at the i-th SF, is 
the sidereal rotational period, and is the number of rotations between Tq 
and Ti. 9q and 9i are the projected directions of the phase angle bisector 
{PAB) in the plane that is perpendicular to the direction of total rotational 
angular momentum at Tq and Tj, respectively. Table 4 shows the epoch of 
SFs. The left hand of Eq. (3) is estimated from the observations; the right 
hand is theoretically calculated from the tentative direction of total rotational 
angular momentum and the orbital information of 107P. We define 5 with 
the following equation 



where is the number of the epoch; here N = A. We can estimate the 
direction of total rotational angular momentum by seeking the minimum of 5. 
Fig. 4 shows the 5 map that is obtained by scanning the celestial sphere with 
a trial axis in steps of 1° in ecliptic longitude and latitude. Two candidates 
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Fig. 4: 5 maps for the direction of total rotational angular momentum by the epoch 
method. 



are found near the directions, A (A = 310°, (3 = -10°) and B (A = 132°, 
/3 = —17°). Since the epoch method generally derives two solutions with 
around the same significance level, we cannot determine a unique solution. 



Table 


4: Epoch of Standard Feature (SF) 


and the 


amount of phase- shift. 






Year/Mon/Day 


Ecliptic longitude 


Ecliptic latitude Amount of phase-shift 




(PAB)[°] 


(PAB)[°] 


(xlO-2) 


2009/11/5 


341.958 


4.611 





2009/12/10 


30.902 


2.805 


-5.305 


2009/12/20 


54.706 


1.094 


-6.113 


2010/01/3 


57.829 


0.879 


-7.794 


2010/02/5 


76.114 


-0.202 


-8.506 



The lightcurve inversion method derives the most adequate shape model 
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Fig. 5: Deviation maps for the direction of total rotational angular momentum by the 
lightcurve inversion method. 



and the corresponding direction of total rotational angular momentum. 
When carrying out the lightcurve inversion method, we set the initial 
conditions for the direction of total rotational angular momentum and the 
sidereal rotation period to that of 0.2979 day. We seek the least deviation 
between the observational lightcurve and the reconstructed lightcurve from 
the shape model by scanning the direction of total rotational angular 
momentum in steps of 1° i n ecliptic longitude and latitude. We use the 
data of November 16, 1979 (IHarris and Young] . Il983l ): December 1-2, 1992 
(jOsip et al.l . Il995l ): and our high photometric precision data (the error is less 
than 0.05 mag), which was obtained from November 7, 2009 to January 18, 
2010. Fig. 5 shows the deviation map of the direction of total rotational 
angular momentum using the lightcurve inversion method. In addition to 
the candidates of the epoch method, we have found three other candidates: 
C (A = 330°, 13 = -27°), D (A = 328°, (3 = -61°), and E (A = 167°, 
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P = 7°). However, we exclude the candidates D and E because they are 
less compatible with the epoch method. Note that the observational data 
was obtained in the phase angle from 21° to 77°. There is no data in a low 
phase angle. Furthermore, the direction of rotational angular momentum 
has an uncertainty of typically more than 5°. For example, though the 
ground-based observati on of Itokawa showed th at the pole orientation was 



A = 355°, (3 = -84° flKaasalainen et al. 



20031 ). the Hayabusa spacecraft 



revealed that the po. 



(jPemura et al. 



e orientation of Itokawa was A = 128.5°, (3 = —89.66° 
20061 ). Moreover, we can see from Table 4 that the SF in 
the lightcurve shifts to the negative direction with time. As a corollary, the 
direction of total rotational angular momentum of the three candidates is 
south of the ecliptic plane. This indicates that the sense of sidereal rotation 
is retrograde. The lightcurve of Fig. 6 has been calibrated for the phase-shift. 
The six peaks and the periodicity in the lightcurve become clearer than in 
the lightcurve of 0.2979 day period in Fig. 2. This result adds to the evidence 
of the retrograde rotation. 

Next, we make the three shape models of 107P for the directions of total 
rotational angular momentum A, B and C. The shape model A is shown 
in Fig. 7. Here, Li, L2, and L3 are defined as the normalized axis length 
when 107P is a triaxial ellipsoid body. The axes satisfy the relationship Li 
< L2 < L3. L3 is a sidereal rotation axis of the shape model A. We have 
found that the normalized axis lengths Li, L2, and L3 are around 1.0, 1.0, 
and 1.6, respectively. The axis ratio of the shape model B is around the 
same as that of the shape model A. The shape model A and B indicate a 
so-called long axis mode (LAM). Some previous studies have described the 
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Fig. 6: Calibrated lightcurve for phase-shift. The fitting curve is described by two-order 
two-dimensional Fourier series. 




Fig. 7: Shape model A of 107P. (Left) Pole-on view. (Center) Equatorial view from the 
right side of the pole-on image. (Right) Equatorial view from the bottom of the pole-on 
image. 
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motion of a force- f ree as ymmetric rigid body ( ISamarasinha and A'Hearru . 
199ll : iKaasalaineru . 1200 ll ). Now, Li ~ L2 indicates that the equations of 
force-free precession are simphfied to the following, 




(5) 



M 



(6) 



/i' 




(7) 



= ii 
20 



{Ll + Ll). 



(8) 



Here, ijj, (f), 6 are the Euler angles of sidereal rotation, precession, and 
nutation, respectively. M is the total rotational angular momentum in an 
inertial frame. Ii and /s express the inertia moment of a triaxial ellipsoid by 
using mass fi. Moreover, the equations show that the motion of the external 



the nutation is negligible, and the angle 6 is constant around 65°. A tilted, 
rugby-ball-shaped body rotates with a period of 0.0993 day about the total 
rotational angular momentum, and with a period of 0.2979 day about the 
external axes of 107P itself. Alternatively, we can assume a case that has 
the sidereal rotation of 0.0993 day and the precession period of 0.2979 day. 
Substituting = 3(f), there is no solution for the nutation angle. Therefore, 
the assumption is not adequate. 

Meanwhile, as we show in Fig. 8, the normalized axis lengths for the 
shape model C are around 1.5, 1.5, and 1.0. Here, the axes satisfy the 



axis about M occurs as a constant rate. From these equations and </> = 31^, 
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relationship Li > L2 > L^. L3 is a sidereal rotation axis. Thus, the shape 
model C is a short axis mode (SAM). The rotation and precession of a SAM 
are in opposite direction from Ji < I3. We calculate the nutation angle by 
assuming <p = —Stp. However, there is no solution for the nutation angle. To 
obtain the solution for the nutation angle, (p should be less than —3.62ip if 
the axis lengths of the shape model C are correct, or the axis lengths of Li 
and L2 should be longer than ^/3 L3 if = —Stp and Li ~ L2 are correct. 
These situations are inconsistent with our results. Therefore, 107P of the 
shape model C is a non-precession object rather than a precessional object. 




Fig. 8: Shape model C of 107P. (Left) Pole-on view. (Center) Equatorial view from the 
right side of the pole-on image. (Right) Equatorial view from the bottom of the pole-on 
image. 



3.3. Taxonomic class and rotational color variations 

The taxonomic class and rotational color variations for 107P were 
investigated by a color-color diagram. We note that the classification of 
subclasses, such as B, F, or G-type, is difficult using multiband photometry. 
We conducted the multiband photometry eight times [Phase-l to -8). The 
obtained color-color diagram and the color index are shown in Fig. 9 
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and Table 5. We utilized the z' images of Phase-Q for those of Phase-7 
due to the poor weather conditions in Phase-7. The color-color diagram 
macroscopically shows that 107P is a C-type (including B, F, G-type) 
asteroid. The colors of Phase-2 and Phase-6 indicate typical C-type features 
in the three color indices. The others are slightly reddish features like an 
X-type asteroid in the color index of r' — i' . Though only the g' — r' of 
Phase-?) barely exceed the one-sigma of mean color index in Table 5, it is 
difficult to assert the detection of the rotational color variation due to the 
photometric error. In addition to it, the long observation term of ~0.15 in 
phase (~1.0 hour) for each sequence obscures the detection of rotational color 
variation. In order to confirm the color variations, follow-up observations 
and/or exploration by spacecraft are needed. 
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Tabic 5: Color index of 107P. The observation term (Obs term) of each 
sequence is expressed as the rotational phase in the lightcurve. Since the data 
of the OAO are obtained with three bands, there is no color index of i' — z' . 
Mean shows the arithmetic average and standard deviation of each color index. 





Obs term [Phase] 


Observatory 


gl _ r' 


r' -i' 


i' - z' 


Phase- 1 


0.0365-0.1221 


OAO 


0.409 ± 0.036 


0.217 ± 0.036 




Phase-2 


0.0048-0.1590 


BSGC 


0.462 ± 0.055 


0.141 ± 0.043 


0.039 ± 0.044 


Phase-3 


0.1348-0.2177 


OAO 


0.522 ± 0.037 


0.184 ± 0.036 




Phasc-4 


0.1941-0.3335 


BSGC 


0.451 ± 0.056 


0.190 ± 0.034 


0.018 ± 0.049 


Phase-5 


0.2345-0.3177 


OAO 


0.364 ± 0.040 


0.177 ± 0.039 




Phase-6 


0.3461-0.5045 


BSGC 


0.444 ± 0.057 


0.138 ± 0.045 


0.050 ± 0.043 


Phaso-7 


0.4662-0.5426 


BSGC 


0.382 ± 0.100 


0.173 ± 0.054 


0.021 ± 0.053 


Phase-8 


0.9590-1.0239 


OAO 


0.423 ± 0.040 


0.176 ± 0.040 




Mean 






0.432 ± 0.050 


0.175 ± 0.026 


0.032 ± 0.015 
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Fig. 9: Color-color diagram of 107P. Let ters in the figure re present the taxonomic classes 
of asteroids on the color-color diagram (jivezic at al.L 120011 ). X-type asteroids include E, 
M, and P-type asteroids. The case of low-albedo asteroids indicates P-type. 
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4. Discussion 



4.1. Tumbling 

We discuss the possibility of tumbling. There were some reports 
i n which asteroid li ghtcurves indicated tumbling, e . g., (2 53) Mathilde 



flMottola et al 



buta tis f ISpencer et al 



(3288) £ 




1995 





Harris et al. 



Kryszczyhska et al. 



19991), and (4179) 



19991 ). 



Pravec et al 



(I2OO5I ) assessed the validity of tumbling for these asteroids based on 



whether the lightcurves could be approximated with two-dimensional Fourier 
series and the physical model of tumbling could be constructed. The 
two-dimensional Fourier series is described in the following form 



m r 

Co 

m m 



E 



^ 27ij . 2nj 

Cjo cos — — t + bjo sm -—t 



+ EE 

fc=l j=—m. 

+ Sjk sin 



Cjk cos 



27rj 2TTk 



2'jTj 2'Kk 



Pik 



Pa, 



(9) 



where m is the order, Co is the mean reduced light flux, Cjk and Sjk are 
the Fourier coefficients for the linear combination of the two frequencies P^^ 
and P^^, respectively, and t is the time. Substituting m = 2, P^ = 0.2979 
day, and P^ = 0.0993 day for 107P, we obtain a fitting curve, as shown in 
Fig. 6. The fitting curve adequately reconstructs the trend of the lightcurve. 
However, since the sidereal rotation period and the processing period have a 
commensurability of 3:1, the lightcurve can also be reconstructed using the 
one-dimensional Fourier series of sixth order. As we mentioned in Section 
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3.2, the physical model is possibly constructed using a LAM of Li-.L^.L^^ 
1.0:1.0:1.6, (A = 310°, (3 = -10°) or (A = 13 2°, = -17° ), 9 = 65°, 



Pravec et al. 



( 120051 ) mentions 



0.2979 day, and = 0.0993 day. Although 
that a tumbling asteroid generally does not return to the same orientation 
in any single period, the approximately equal length of Li and L2 for 107P 
suggests a negligible change for the nutation angle. Therefore, 107P can 
return to the same orientation every 0.2979 day. These circumstances imply 
that 107P might be a tumbling object. 

Assuming 107P is a tumbler, external forces are required to trigger the 



motion. Impacts of small objects, tida. 



effects are suggested by 



Pravec et al. 



e ncoun ters with planets, and YORP 
(I2OO5I). Though 107P is a NEO, 



the object did not have an encounter with Earth in 1949. In the case of 
km-size objects, the efficient onset of tumbling b y YORP requires a longe r 
timescale than that of collision with small objects ( IVokrouhlicky et al.l , 120071 ). 
Therefore, we propose the impact of small objects as a probable cause for 
tumbling of 107P. The orbital origin of 107P has a high possibility of being 
from the outer MBA region inhabited by MBCs. One possibility is that the 
cometary activities of MBCs are caused by impacts of small objects. We can 
consider that 107P is originally an object like an MBC and impacts with 
small objects in the NEO region could eject dust and/or expose sub-surface 
ice that then trigger 107P's cometary activity. When we sup pose that the 
collisional excitation happened in 1949, the damping timescale (IHarrisI , Il994j ) 
is expressed as 



p3 



(10) 
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where D is the mean diameter of tumblers in kilometer units and C is a 
constant of about 17 (uncertain by about a factor of 2.5). The units of 
and r are hours and billion (10^) of years, respectively. Since the damping 
timescale of around 6.2 x 10^ yr is long enough, 107P would continue tumbling 
even if the impact occurred before 1949. 

4-2. Binary asteroids 

We describe the situation in which 107P hosts a binary. In order to 
confirm the existence of a binary, the detection of mutual eclipse events is 
required in the lig htcurve. The mutual ecli pse e vents were not de tected in 

fll995l ) because 



the observations of 



Harris and Young! (119831 ) and 



Osip et al. 



of the viewing angle, the lower photometric precision, or the absence of the 
binary. On the other hand, we detected the around same flux decrease in 
every 0.50 phase. Therefore, the existence of the binary is conceivable as 
the other interpretation of the shape model C. If we define the flux decrease 
around the phase of 0.15 (or 0.30, 0.45) and 0.65 (or 0.80, 0.95) in Fig. 
6 as the primary (secondary) eclipse and the secondary (primary) eclipse, 
respectively, the orbital period of the binary is 0.2979 day. Supposing a 
circular orbit and negligible mass for the binary, the semi-major axis is 
described as 



47r2 



where G is the gravitational constant, M is the mass of 107P, and Porb is 
the orbital period of the binary. For the sake of sim plicity, when assuming 
that 107P is spherical with the diameter of 3.46 km (ILicandro et al.l . 120091 ) 
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and a typical density of 2 g/cm?, the semi-major axis is around 3.65 km. 
In the case of the same albedo for 107P and the binary, th e flux decrease o: 



the t otal eclipse {A^ut) satisfies the following relationship f lPolishook et al 



201l|) 



Amut = 2.5log 



Up 



(12) 



where Rg is the radius of the binary and Rp is the radius of 107P. Since 
the typical flux decrease is ~0.05 mag in Fig. 6, the radius of the binary is 
around 0.4 km. When we assume the orbital plane of 107P accords with the 
line of sight from an observer, the inclination of the binary as an occulter 
satisfies 

sim<^^^±^. (13) 
a 

Here, i is the inclination of the binary for the orbital plane of 107P is less than 
36° in the 107P system. If i is zero, the eclipse duration is estimated to be 
~0.05 day. The term is around one-sixth of the orbital period and consistent 
with the interval of lightcurve peaks of Fig. 6. Moreover, the binary 
hypothesis indicates that the double-peak period of the lightcurve without 
the eclipse becomes 0.1490 day. As we mentioned in Section 3.1, the period 
of 0.1490 day as the sidereal rotation of 107P is not likely. Alternatively, the 
lightcurve without the eclipse might be a quadruple-peak lightcurve whose 
period is 0.2979 day. Though the quadruple-peak lightcurve is rare, the 
period could compatibly account for all the past reports. In addition, the 
situation shows that the sidereal rotation of 107P and the orbital periods of 
the binary are locked with 0.2979 day. The period of 0.0993 day is explained 
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by the period between the egress time of the primary (secondary) echpse and 
the ingress time of the secondary (primary) echpse. 

The promising mechanisms for formation of asteroid binaries are the 
rotational-fission due to the spin-up by Y ORP effects (IScheeres 



2007 



Pravec and Harris . 



(IRichardson et al 



2007 



1998 



Walsh et ah 



Bottke et al. 



2008 



1999 



) , tidal encounter with planet s 



Walsh and Richardsonl. 



and t he escaping ejecta by the collisions (IDurda et al. 



2004 



20061) 



Polishook et al. 



201ll ). The mechanisms have a lot in common with the cause of tumbling. 



Fissions and collisions in every mechanism can trigger 107P's cometary 
activity. The possible existence of a binary is consistent with the past 
cometary activity. 



5. Summary 

This study revealed the physical properties of 107P by a photometric 
observation campaign. We detected the lightcurve periodicity to be 0.2979 
day and 0.0993 day with a commensurability of 3:1. The multiband 
photometry indicates that the taxonomy class of 107P is C-type. No clear 
rotational color variations are confirmed on the surface. We suggested two 
models to explain the different interpretations of the lightcurve periodicity. 

1. The commensurability reflects tumbling with the sidereal rotation 
period of 0.2979 day and the precession period of 0.0993 day. The shape 
is a LAM of Li:L2://3 = 1.0:1.0:1.6. Around the same length of Li and L2 
shows the nutation angle is approximately constant at 65°. The direction 
of total rotational angular momentum is around A = 310°, /3 = —10°, or 
A = 132°, fi = —17°. 107P returns to the same orientation every 0.2979 
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day by retrograde motion. Impacts of small objects are suggested as a 
cause for the tumbling and comet activity. Alternatively, the past comet 
activity itself is thought to be a cause of the tumbling, like a IP/Halley 
(ISamarasinha and A'Hearnl . Il99ll ). 

2. 107P is not a tumbler. The sidereal rotation period is 0.2979 day. The 
shape is roughly spherical but slightly hexagonal with a SAM of Li:L2:L^ = 
1.5:1.0:1.0. The pole orientation is around A = 330°, /3 = —27°. The sense 
of rotation is retrograde. The lightcurve of commensurability would reflect 
a discriminative appear ance like (28 6 7) Ste ins, which has been explored by 



the Rosetta spacecraft (IKeller et al 



20101 ). Otherwise, the lightcurve also 



indicates the possibility of hosting a binary whose orbital period is 0.2979 
day. The existence of a binary is also consistent with the past cometary 
activity. 

Finally, we describe the mission feasibility for 107P. The orbit accessible 
by spacecraft makes 107P a promising target for a sample-return mission. 
If 107P is not a tumbling object, the moderate rotational period of 0.2979 
day would enable us to obtain a sample by the touchdown of a spacecraft, 
whereas touchdown on 107P would require a difficult maneuver if 107P is a 
tumbler. In that case, a multi-fiy-by mission that combines with the sample 
return mission for another target would become a hopeful plan. 
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